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ABSTRACT 

We study the far-infrared emission properties of the nearby starburst galaxy 
NGC 253 based on IRAS maps and an ISOPHOT map at 180 fim. Based on the 
analysis of the light profiles, we have been able to identify three main structural com- 
ponents: an unresolved nuclear component, an exponential disk, and a kiloparsec scale 
bar. In addition, we also found a ring structure at the end of the bar that is particularly 
conspicuous at 12/xm. The Spectral Energy Distribution (SED) of each morphological 
component has been modeled as thermal dust emission at different temperatures. The 
unresolved nuclear component is dominated by cold dust emission (T,„ ~ 50 K), whereas 
the disk emission is dominated by very cold dust (T^^ ~ 16 K) plus a contribution from 
cold dust (T^ ~ 55 K). The bar emission corresponds mainly to cold dust (T^ ~ 23 K) 
plus a warm component {T'^ ~ 148 K). We detect an extension of the disk emission due 
to very cold dust, which contributes a large fraction (94%) of the total dust mass of the 
galaxy. The estimated total dust mass is 8.2 x 10'' Mq. 

Subject headings: galaxies: individual (NGC 253) - galaxies: starburst - galaxies: 
photometry - galaxies: structure - infrared: galaxies - interstellar dust. 



1. Introduction 



Much of our knowledge about the infrared (IR) emission from galaxies comes from their total 
emission fluxes. Models of the spectral energy distribution (SED), from UV to radio wavelengths 
have been successful in reproducing the observations in a variety of galaxies, from normal to star- 
burst types (Devriendt, Guiderdoni, & Sadat 1999; Silva et al. 1998; Efstathiou, Rowan-Robinson, 
&: Siebenmorgen 2000). Nevertheless, there are still some controversial aspects which need to be 
addressed, such as: 

Whether heating of dust by ultraviolet photons from young stars is the dominant source of 
the far-IR (FIR) emission (Devereux et al. 1994), or if there is a noticeable contribution 
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from non-ionizing photons from bulge and disk stars (Walterbos & Schwering 1987; Xu & 
Helou 1996; Mayya & Regaradjan 1997). 

The inner gas-to-warm dust mass ratios in external galaxies result in about 5-6 times the 
local Galactic value (Dcvcrcux & Young 1990), implying a deficit of cold dust. This can 
be attributed to a lack of data beyond 100 fiui, yielding both an overestimate of the dust 
temperature and an underestimate of cold dust masses hidden by other, more energetically 
important, components. 

The extension of dust disks as compared to stellar disks observed in the optical bands. En- 
gargiola (1991) and Engargiola & Harper (1992) found that the scalelength of stellar disk 
is similar or larger than the observed at 160/Lim continuum. However recent work seems to 
support dust disks larger than the stellar ones by factors from 1.4 up to 2 (Davies et al. 1997; 
Alton et al. 1998; Xilouris et al. 1999). 

What is the opacity of the disk in spiral galaxies? Disney, Davies & Phillipps (1989) sug- 
gested that the spiral galaxies could be optically thick. However, more recent works indicate 
only moderately opaque disks (Xilouris et al. 1999). This feature has relevant cosmolog- 
ical implications, giving different degrees of extinction to be accounted for in the cosmic 
background. 

Adding spatial resolution to the study of the far-IR emission from galaxies certainly helps re- 
solving most of these controversial points. This has been possible only in a limited number of nearby 
galaxies thanks to the instruments on board the NASA-Kuiper Airborne Observatory (KAO) and 
to the improved resolution ('^ 1 arcmin) of maps obtained with the Infrared Astronomical Satellite 
(IRAS) and processed by the HiRes algorithm (Aumann, Fowler, & Mclnyk 1990; Rice 1993). 
Engargiola (1991) has studied the FIR continuum distributions of NGC 6946 from 60 to 200 /im, 
using IRAS and NASA-Kuiper Airborne Observatory data, showing that the light distribution re- 
veals discrete sources coinciding with giant H II regions superimposed on an unresolved exponential 
dust emission background. The brightness profiles were separated into a nuclear starburst plus 
an exponential disk. Both the nuclear and the disk SEDs were described as modified blackbody 
sources at T ~ 40 K and T '-^ 15 K, respectively. Far IR maps of the quiescent spiral galaxies 
M 31 and NGC 4565 have been studied by Xu k Helou (1996) and Engargiola & Harper (1992), 
respectively. In both cases they could distinguish warm dust tracing the spiral arms superimposed 
on a diffuse cool dust emission corresponding to the disk. They also found that old stars contribute 
an important fraction of the heating of the dust. 

More recently, with the launch of ISO with its on-board instrument ISOPHOT, maps of similar 
or better resolution have been obtained up to 200 /xm, allowing more detailed studies of various 
physical quantities related to the FIR emission of galaxies. Questions concerning the dust and tem- 
perature distributions, the structural components followed by the dust distribution, the existence 
of very cold dust emission, and the physical extent of the disks can all be addressed. Haas et 
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al. (1999) have mapped the FIR emission of M 31 at 175 /xm using ISOPHOT data, showing the 
importance of the contribution of cold dust to the total mass of the galaxy. 

Relevant advances also come from the theoretical front. Most current models are able to 
reproduce the main features of the far IR emission from different types of galaxies consistently with 
other wavelength ranges dominated by the stellar radiation. For instance, Devriendt, Guidcrdoni, 
& Sadat (1999) developed successful models which take into account the evolution of the stellar 
population as well as extinction and emission of dust in a self consistent way. On the other hand, 
Bianchi, Davies, k, Alton (2000) constructed bidimensional models to predict the dust-reprocessed 
FIR output in spiral galaxies corresponding to different stellar and dust geometries. These models 
could explain the observed far IR spectral energy distribution from spiral galaxies only when dust 
disks are optically thick. However they failed to reproduce the large far IR scale lengths suggested 
by recent far IR observations (Alton et al. 1998). 

Because of its proximity, NGC 253 is a good candidate for detailed studies at any wavelength 
range. Its angular size in the visible range is 27.5' x 6.8'. It is classified as a spiral galaxy of 
type SAB(s)c and is seen almost edge-on. Throughout this paper, we adopt a distance of £) = 3.4 
Mpc (Sandage & Tammann 1975), implying a scale of 16.5 pc/arcsec. NGC 253 is considered as 
a prototype nuclear starburst galaxy, exhibiting strong IR emission {Lir ~ 2 x lO^^Lo), being an 
ideal target for IR studies. In this paper, we analyze maps of NGC 253 from 12 jiva. to 180 iim.. 
We are able to separate the far IR emission into its morphological components. The SED of each 
component was been modeled in terms of blackbodies emission. 

2. Mid and Far IR maps. 

The ISOPHOT (Lemke et al. 1996) data at 180 //m were retrieved from the ISO^ Data 
Archive. The map was performed in undersampled raster mode with a step of 180 arcsec, using 
the C200 camera, which consists of a 2 x 2 pixel array (the pixel size is 89 arcsec). The data 
were reduced using PIA v9.1 (Gabriel et al. 1997), and the following corrections were applied: 
ramp linearization, deglitching at the ramp and signal levels, and dark signal subtraction. Only the 
last 60% of the signal at each raster position was retained in order to minimize transient effects. 
The flux calibration was done using the internal calibration lamps, although the final calibration 
was determined by comparing the sky background measurement in the 180 /xm ISOPHOT map 
with background measurements from C05£'/DIRBE maps. A scaling factor of 1.48 has to be 
applied to the ISOPHOT data in order to match both numbers. Two pixels of the C200 array 
were saturated when crossing the nucleus of the galaxy. The values at these map positions were 
obtained by interpolation from the neighbourings . This fact may produce a broadening of the 



^Based on observations with ISO, an ESA project with instruments funded by ESA Member States (especially the 
PI countries: Prance, Germany, the Netherland and the United Kingdom) and with the participation of ISAS and 
NASA. 
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light profile at the center of the galaxy, and probably underestimate the nuclear contribution. The 

map reconstruction was performed within PIA using the trigrid algorithm. We estimate the actual 
resolution of this map to be ^ 180 arcscc, due to the poor raster sampling. The isocontours of the 
180 lira image overlaid on a i? band (MPG/ESO 2.2 m + WFI optical image^) are presented in 
Fig. 1. The main features present in the optical image can also be recognized in the 180 ;um map. 
However, the morphology at 180/xm is slightly different than that observed in the other IR bands. 
The emission along the minor axis extends Tip to 10 Kpc from the nucleus. This extended emission 
can be associated with a cold dust halo. Addittional evidence was found in our analysis of the strip 
maps perpendicular to the disk (Perez Garcia et al. 2001) The same conclussion has been reached 
by Radovich, Kahanpaa, & Lemke (2001). A similar morphology has been reported from ROSAT 
X-ray observations by Pietsch et al. (2000). 

We retrieved IRAS maps at 12, 25, 60, and 100//m from the Infrared Processing and Analysis 
Center at Caltech (IPAC). The data retrieved were processed using HiRes algorithm, which is based 
on the Maximum Correlation Method (Aumann, Fowler, &; Melnyk 1990) to produce images with 
improved spatial resolution. The resolution of these IRAS maps ranges from 1 arcmin at 12 and 
25 fiui to 1.7 arcmin at WOfim. Unfortunately, the low level emission seen in the IRAS maps 
cannot be used due to the spurious structures caused at 12 and 25 fim (see Fig 2) by detector 
hysteresis after crossing the bright nucleus, and at 60 and 100 /im by reflected emission off the 
IRAS secondary spider (Rice 1993). 

Very recently, Radovich et al. (2001) have studied the far IR emission of NGC 253 using IRAS 
and ISO data, as well as a numerical radiative transfer model. Their main aim war detecting FIR 
emission above the disk that might be connected with outflows and/or supergalactic winds. These 
authors have also shown the need for an extended dust disk. 

3. Results 

3.1. IR luminosity profiles 

For all maps, we obtained surface brightness profiles by azimuthally averaging over elliptical 
annuli. We fitted the isocontours of surface brightness with the ELLIPSE task in IRAF^, which 
uses the algorithm described in Jedrzejewski (1987). The spatial interval between two consecutive 

isophotcs is equal to the pixel size in each filter, but averaged over an annulus of width close to 
the resolution in each filter. Figure 2 shows the family of ellipses fitted for two filters, 12 and 25 
/xm. At first glance, it can be noticed how the isophotal twisting delineates the presence of a bar 



^Obtained with the WFI (full field) at the MPG/ESO 2.2 m telescope and retrieved from the ESO archive. 

^IRAF is the Image Reduction and Analysis Facility, a general purpose software system for the reduction and 
analysis of scientific data. TRAP is written and supported by the IRAF programming group at the National Optical 
Astronomy Observatories (NOAO). 
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and a ring at the circumnuclear region (discussed below). The brightness profiles, ellipticities (e), 
and position angles (PA) in all filters as a function of radial distance from the galaxy center are 
shown in Fig. 3. Fig. 4 shows the light profiles normalized to the central peak to facilitate the 
comparison at the different wavelengths. The brightness profiles show similar behavior in all the 
IR filters, with the exception of the 180 //m profile. Two different slopes are clearly distinguished, 
a steep profile from the central peak out to about 2-4 arcmin, and from there on the profile flattens 
sampling the disk of the galaxy (see Fig. 4). The IR emission concentrates at different radii 
depending on the wavelength: the sharpness of the luminosity profiles from 12 to 100 fim reflects 
the predominance of the warm dust emission associated with the starburst activity at the galaxy 
center. In the outer parts the disk is the dominant structure. In contrast, the luminosity profile at 
180 /xm appears flatter than the other profiles, following an exponential disk profile behavior at all 
radii. No central region can be distinguished in the 180/im profile. The absence of a central peak 
at 180//m is probably a consequence of the lower prominence of the nucleus at that wavelength, 
although the saturation problems mentioned above likely blur the central peak. Remarkably, the 
extension and the decay of the light profile at 180 /im is comparable to that measured in the B 
band (Pence 1980) (see Fig 4), although the slope beyond 15 kpc is slightly flatter at 180 /xm than 
at the B band. This fact may indicate an extended disk of very cold dust (Radovich et al. 2001). 

More structures can be distinguished in the brightness profiles. At a distance of ~ 3.5 kpc, 
there is a shoulder in the 12 /xm profile, which is also evident in the 25 /xm profile; this feature 
coincides with the end of the inner disk detected by Scoville et al. (1985), which we have also 
identified as a ring structure (discussed in section 3.2). There are other findings which support 
this identification: an excess of millimeter CO emission at 3.5 kpc along the semi-major axis, 
reported by Scoville et al. (1985), and more recently by Sorai et al (2000); and a Ha image 
shows a remarkable ring structure with radius of 2.7 arcmin (Hoopes et al. 1996). The prominent 
hump observed in the B band at a distance between 6 and 12 arcmin along the semi-major axis is 
marginally present in the 100 /xm profile. The hump at the B band has been attributed to young 
stars and Hll regions in the spiral arms (Pence 1980). Note that in the other IRAS filters we cannot 
see the hump, probably because of noise introduced by spurious structures at low flux levels. 

The position angle and the ellipticity also vary considerably with galactoccntric distance. In 
particular, the 12 iim light profile shows drastic changes in both ellipticity and position angle at 
radii between 2 and 3.5 kpc. The ellipticity increases from the center up to a maximum value of 
0.8 at a distance of 3.5 kpc, where it suddenly decreases to a value of 0.5; from there outwards, the 
ellipticity approaches the disk value of 0.75. Simultaneously, the position angle changes from 70° 
to 60° at 2 kpc falling to the outer disk average value of 52.3° at a distance of 3.5 kpc. The 25 /xm 
profile show a similar behavior to that at 12 /xm, although the position angle remains constant at 
70° up to 3.5 kpc. We identify these features with the presence of a bar, whose existence has been 
reported previously by other authors (Pence 1980; Scoville et al. 1985; Forbes & DePoy 1992). 
Both the orientation and the extension match well the values reported in the near IR: Scoville et al. 
(1985) found in a 2.2 /xm map an elongated barlike structure with a radius ~ 1.98. kpc, oriented 



-6- 



at PA = 68°, within an inner disk which extends out to ~ 2.97 kpc; Forbes &; DePoy (1992) found 

in an H band image a bar oriented at ~ 70° and extending 2.48 kpc. The barred nature of this 
galaxy is obvious in the near IR (Scoville et al. 1985) and also quite evident in the 2MASS Atlas 
Image Gallery mosaic of this galaxy^. Moving outwards, the position angle and ellipticity reach 
the same asymptotic limiting values for the different filters, within the errors. We obtain for the 
ellipticity the average value 0.79 ±. 0.10, which corresponds to a galaxy inclination of i ~ 78° it 6°, 
and position angle 50° it 7, representing the disk. These values are in good agreement with those 
obtained by Pence (1980) using optical data (PA = 51°, i ~ 78.4° ± 0.3°). As the resolution 
decreases — moving to longer wavelengths — details of the isophotes at the inner part become less 
representative. At 180 /xm, the position angle reaches similar values in the outer parts to those at 
12, 25, 60, and 100 jim.. However the ellipticity remains lower, the isophotes being rounder. This 
is consistent with the detection of a cold dust halo. This halo is not detected at other wavelengths 
except at 100 //m which shows a slight lower ellipticity. This is consistent with the work by Alton 
et al. (1998) who fail to detect dust emission above the disc using IRAS images from edge-on 
galaxies. 

We have also looked at the variation of the IR colors with galactocentric distance (see Fig. 
5). The mid- and far-infrared color ratios approximately describe the origin of the dust emission 
at these wavelengths (Helou 2000; Telesco 1993). In a simplified scenario, the IR emission of a 
relatively normal galaxy can be described as a composite of two spectral components (Helou 1986): 
blackbody emission from classical grains in temperature equilibrium heated by the interstellar 
radiation field dominating in the far IR; and very small grains and PAH molecules transiently 
heated to high temperatures by the absorption of single energetic photons produced in star forming 
regions. As a result starbursts have warmer 60/100 /^m color temperatures and cooler 12/25 ^m. 
color temperatures. The more active starbursts exhibit bluer 25/60 /xm and 60/100 fim. colors 
corresponding to warmer dust in high temperature equilibrium. As a guide, we have superposed 
for all colors the typical values of two model components — the disk (D) and the starburst (SB). 
These models were used to synthesize the observed far-infrared spectra of many IRAS galaxies by 
Rowan-Robinson & Crawford (1989). The first thing to notice is that in the outer regions all 
colors are typical of disk emission. The 12/25 jJLm color index shows a strong peak, coincident with 
well identified morphological features (described in next section). In the innermost regions, this 
color matches that of typical starburst galaxies. Around 3-4 kpc there is a strong color increase 
slightly beyond the end of the bar at the location of the ring detected in our profile separation 
(as described below). The origin of this color enhancement can be related to an intense PAH 
emission associated with the ring of molecular gas mentioned above (Scoville et al. 1985; Sorai 
et al 2000). The 25/60 //m color index shows a smoother trend. An interesting feature appears 
around 2 kpc, coincident with the head of the barlike structure, that might be another indication 
that star formation is taking place there. The 60/100 /xm and 100/180 /xm color indexes have an 



''This image (J,H & Ks composite) is available electronically at: http://www.ipac.caltech.edu/2mass/gallery/- 
images_galaxies.html. 
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absolute maximum at the center of the galaxy, as corresponds to the central starburst activity. 
The 60/100 nm color index shows a small increase coincident with the beginning of the spiral arm. 
This is interpreted as an increase in dust temperature due to star formation in large HII regions, 
as observed by Engargiola (1991) in NGC 6946. 



3.2. Surface brightness profile decomposition 

The best fit structural components to the radial light profiles at 12, 25, 60, 100 &: 180//m are 
shown in F. 6. Initially, we tried a morphological separation of the profiles into two components: an 
unresolved source corresponding to the PSF (for the IRAS maps we use gaussian functions and for 
the ISOPHOT map a model PSF) and an exponential disk that mainly accounts for the external 
regions of the galaxy. The resolution of the 180 fim map is not sharp enough to separate any other 
components than the exponential disk. A bulge component was not included because it cannot be 
sampled with the limited resolution of our maps. Forbes & DePoy (1992) found an effective radius 
(10") for the bulge which is well below our resolution limit, so its contribution will be implicitly 
included in the unresolved component. The analytic profiles were convolved using the corresponding 
PSF function for each filter before performing the fit. A third component was introduced to improve 
the matching of the models with the observed profiles following the indications of the presence of 
a bar discussed in Section 3.1. For this latter component we have taken a fiat bar profile (Prieto 
et al. 1997): 

T / \ -^bar 

-'barv J 7 \") 

l + exp(^) 

where /bar is the amplitude, rtar is the length and ^bar is the downward gradient. 
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Table 1. Structural parameters 



A 


IpSF 


^disc 




^bar 


^bar 


^bar 








/im 


MJy/sr 


MJy/sr 


kpc 


MJy/sr 


kpc 


kpc 


MJy/sr 


kpc 


kpc 


12 


100±20 


18±2 


3.2±0.3 


280±40 


0.73±0.08 


0.26±0.03 


8.5±0.9 


3.0±0.3 


0.69±0.11 


25 


1950±200 


19±2 


3.7±0.4 


190±20 


1.34±0.14 


0.28±0.03 


0.3±0.4 


3.0±0.3 


0.69±0.11 


60 


7880±790 


99±10 


4.3±0.4 


130±20 


2.1±0.2 


0.35±0.12 


0.2±3.9 


3.0±0.3 


0.69±0.11 


100 


2760±280 


230±20 


3.5±0.4 


550±70 


1.9±0.2 


0.37±0.04 


4.8e-4±2.4e-4 


3.0±0.3 


0.69±0.11 


180 




550±170 


3.1±0.9 
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In addition, the 12 fim luminosity profile shows a conspicuous maximum around 3.5 arcmin, 

at the position where a strong decrease in the ellipticity and a noticeable variation in the position 
angle are found. We identify this feature with a ring at the end of the bar, which also coincides 
with the maximum extent of the inner disk identified by Scoville et al. (1985). This structure can 
be related to the molecular and Ha rings detected by(Sorai et al 2000) and (Hoopes et al. 1996). 
A ring structure can be described by a Gaussian profile (Buta 1996), as follows: 



where /ring is the amplitude, rring is the center and Zring is the width. 

The ring can still be recognized in the intensity profile at 25 /xm, but it disappears at 60, 100 
and 180 /im. Nevertheless, wc have introduced this morphological component in the decomposition 
of the 25 and GOfim profile, although the radial distance (rii^g) and the width (/ring) are fixed to that 
at 12//m, in order to obtain a meaningful fit. We have checked that introducing this component 
at the other wavelengths (60 and 100 /xm) the fit results do not change, being compatible with a 
very low amplitude for it. The uncertainty of the amplitude of the ring at 60/xm is compatible with 
non-detection (see Table 1) the seeing ring. 

The resulting parameters for the different morphological components and filters are given in 
Table 1. During the fitting process all parameters were allowed to vary freely, except for the radial 
distance and width of the ring, as mentioned above. The uncertainties in the parameter values were 
computed using a bootstrapping technique. This procedure is based on a Monte Carlo simulation: 
new light profiles are obtained by perturbation of each measTircd profile value, using a normal 
distribution of the same width as the error of that point. New fit parameters are determined for 
each simulated profile (usually 50 simulations) and the uncertainty of each parameter is taken as 
the standard deviation of the resulting values. 

The disk scale lengths obtained in all the examined maps are very similar and have a mean 
value of 3.63 kpc. This is very close to the value of 3.32 kpc obtained in the B band by Pence 

(1980). Puche & Carignan (1991) found in H l observations a disk scale length of 3.41 kpc. 
However, Forbes & DePoy (1992) obtained in the H band a disk scale length of 2.65 kpc along 
the semi-major axis. The smaller scale length observed in the H band, if confirmed at ./ and K, 
could be attributed either to high extinction at the inner parts or to the presence of very high 
temperature small dust grains emitting in the near IR. An important result is that both the stars 
scale length (Pence 1980) and the dust scale length are similar, which is not in agreement to the 
results finding of Radovich et al. (2001) based on radiative transfer models of this galaxy. 
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3.3. Spectral energy distributions of the structural components 

The IR spectral energy distribution (SED) of the morphological components, namely, the 
unresolved nucleus, the disk, and the bar found in the previous section, are displayed in Figure 7. 
For each filter and each morphological component we represent the emission integrated over the 
galaxy. The uncertainty of the emission from any component is determined from the uncertainty 
of the amplitude. The global emission of each morphological component has been modeled as the 
combination of the emission from emissivity-weighted blackbody components (Hildebrand 1983). 
An emissivity index of 2 has been adopted. The best fitting blackbody temperatures are compiled in 
Table 2. The uncertainties in the temperature and scaling factor of each blackbody were computed 
using the bootstrapping technique, similar to what was done for the profile decomposition. 

The SED of the unresolved nuclear component can be modeled as a cold (T ~ 49 K) compo- 
nent. The cold component has a temperature typical of dust heated in star forming regions. This 
component is present in nearly all classes of galaxies, including normal, starburts, and even active 
galaxies (Knapp et al. 1996; Chini, Kruger, & Kreysa 1992; Klaas et al. 1997; Perez Garcia 
& Rodriguez Espinosa 2001). There is an excess at 12^m which is likely due to a combination 
of PAH and very small grains emission. In fact, PAH emission has been detected at the central 
regions of NGC 253 (Perez Garcia et al. 2001; Dudley & Wynn-Wilhams 1999; Keto et al. 1999) 

The SED of the disk component is explained by a combination of a cold (T ~ 55 K) and a 
very cold (T 2± 15 K) component. The origin of the warm component is again due to dust heated 
in H II regions and OB associations populating the arms of the disk. The temperature of the very 
cold component is typical of dust heated by dilute interstellar radiation. This very cold dust has 

been observed in normal galaxies (Walterbos & Greenwalt 1996; Walterbos & Schwering 1987; 
Cox, Kruger, & Mczger 1986) and in the disks of starbursts and active galaxies (Radovich et al 
1999; Perez Garcia & Rodriguez Espinosa 2001). 

The SED of the bar component can be modeled as the sum of two modified blackbodies, a 
warm one at T = 148 K and a cold one at T = 22 K. The origin of the warm component is related 
to regions of very intense star formation Klaas et al. (1997); Lutz et al. (1996). Instead, the 
temperature of the cold component indicates an origin at a quiescent interstellar gas region, where 
star formation activity is low. 

Most of the emission from the ring is concentrated in the mid-IR range, from 12 to 25 /itm; 
beyond this wavelength the emission cannot be reliably determined. The characteristics of the 
ring emission seems to indicate a dominant contribution of PAH and very small grains, as already 
mentioned in Section 3.1, indicative of photodissociation regions. 
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4. Derived physical properties 

4.1. IR luminosities 

We have computed the IR luminosities between 1 and 1000 fim for each morphological compo- 
nent (except for the ring because of the reasons mentioned in the previous section). The luminosity 
for each morphological component can be easily computed after the blackbody model is found. The 
resulting values are compiled in Table 2, the total values are obtained as the sum over all structural 
components. 
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Table 2. Temperatures, infrared luminositites, dust masses and star formation rates derived for 
the different components, as well as the global values for the galaxy. 



Component 


T(K) 


Fir 




Dust mass 


SFR 




(K) 


[10-® erg s-l cm-2] 


[109 Lq] 


[10« Mq] 


(Mo/yr) 


PSF 


49.4±0.9 


5.77 ± 0.03 


20.9 ± 1.5 


0.21 ± 0.02 


3.5±0.3 


DISK 


55±15 


1.097 ± 0.011 


4.0 ± 1.3 


0.02 ± 0.03 


0.7±0.2 




15.9±2.1 


2.1 ± 1.6 


8 ± 3 


77 ± 58 




BAR 


148±74 


0.8 ± 0.5 


3 ± 2 


« 1 






23±12 


1.4 ± 0.9 


5 ± 3 


5 ± 14 


0.8±0.6 


TOTAL GALAXY 




11.4 ± 1.5 


41.1 ± 1.6 


82 ± 7 


5±1 
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The nuclear unresolved component is the most important contributor to the total IR luminosity, 
accounting for about half of the total emission. The disk accounts for ~ 30% of the total, of which 
2/3 corresponds to the very cold component. The remaining ^ 20% comes from the bar structure, 
in which the cold component dominates. The contribution of the ring to the total luminosity is 
very small. 

4.2. Star formation rate 

Star formation rates (SFRs) were calculated following Kennicutt (1998) for both the warm 
and cold components using the following expression: 

SFR = 1.7 X 1O-^°Lir/L0 [Mq yT~^]. 

This relationship applies only to relatively young starbursts, which implies that it only produces 
meaningful results when T^ust ^ 25 K. The resulting values arc compiled in Table 2. As expected, 
the maximum star formation activity is taking place in the nuclear region at a moderate rate 
~ 3.5 Mq/jy. Furthermore, star formation activity at a lower rate is taking place in both the 
disk and the bar. The total SFR in the galaxy is 5.0 MQ/year. This value is 20% lower than 
that found by Radovich et al. (2001) who give a total SFR value of 6.1 Mq/ji (after correcting 
their value for the different distance used herein). Notice that Radovich et al. (2001) include the 
20 K dust emission to complete the global SFR. We therefore consider the value here provided as 
a rather more precise value. If we compare this value (4.97 M0/year) with that obtained from the 
Ha emission, SFR = 0.6192 Mq/jt by Hoopes et al. (1996), we estimate that the Ha emission 
must be obscured by a factor of ~ 8, in good agreement with the value estimated by Hoopes et al. 
(1996). 

4.3. Dust masses 

Dust masses were estimated following Hildebrand et al. (1977) for the emissivity index (/? = 2) 
used here, assuming dust grain properties as in Hildebrand (1983): 

Md = 3.6 X 10-5(ri^/40)-6LiR/LQ [Mq] 

The computed values for each component are shown in Table 2. The total dust mass of the 
galaxy is computed as the sum of all contributions. The most important contribution (about 94%) 
to the total mass comes from the very cold dust present in the disk. The mass estimate depends 
on the emissivity law adopted; we have checked that for the emissivity index 1, the resulting mass 
varies by a factor 0.7 to 2 for dust temperatures in the range 18 to 65 K, respectively. We have 
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also investigated the radial distribution of the dust mass by integrating over eUiptical annuh at a 

given radius. The cumulative distribution per component is plotted in Figure 8. It can be noticed 
that about half of the total dust mass of the galaxy is contained in the first 4-5 arcmin. 

The total dust mass value thus estimated gives a more reliable value than that obtained from 
the global SED. As an exercise, we have computed the mass using the integrated emission of 
the galaxy. For that purpose, the SED have been modeled as dust emission from a cold (T ~ 
54 K) and a very cold (T ~ 20 K) component. In this way, we derive a total luminosity of 
[3.68 ± 0.02] xlO^° Lq and total dust mass of [2.92 ± 0.04] xlO^ Mq, which are in agreement 
with the total values obtained for the whole galaxy by Radovich et al. (2001) (once the difference 
in distance has been allowed for). Note that this value is ^ 3 times smaller that the dust mass 
obtained with the multiple-component approach. 

We have also computed the gas-to-dust mass ratio for diff'erent galactocentric distances using 
gas mass data from the literature (Scoville et al. 1985). The result is plotted in Figure 9. The 
mass ratio decreases strongly towards the outer part of the galaxy approaching the canonical value 
of our Galaxy. In our Galaxy the gas-to-dust mass ratio is ^ 160 (Sodroski et al. 1994). Devereux 
k, Young (1990) obtained a value of 1070 for the inner disk (the inner half of the optical disk, 
6.8 kpc) of NGC 253 using IRAS data. They concluded that this value is too high and would 
indicate that ~ 80-90% of the dust mass in spiral galaxies is radiating at A > 100 /xm. At 5.2 kpc 
from the center, using our dust mass estimate, we obtain a value of ^ 60. Furthermore, Houghton et 
al. (1997) provided values for molecular and atomic hydrogen for R < 5.6 kpc, the total hydrogen 
mass being 2.4 x 10^ Mq. We find that the global gas-to-dust mass ratio is ~ 30, which is close to 
the outermost value in Figure 9. 

Summarizing, we find that the gas-to-dust mass ratio varies greatly along the galaxy, from 
a value of 300 at 1 kpc to ~ 60 in the outermost part. The dust mass appears to be distributed 
towards the outer parts of the galaxy, and is rather difficult to detect at A^lOO/xm due to the very 
low temperature. The gas to dust mass ratio seems to indicate little excess of dust as compared to 
the canonical value in our Galaxy. 

5. Summary and conclusions 

We have used ISOPHOT at 180/im and IRAS maps to study the spatial distribution of the 
far-infrared emission in NGC 253. We have performed an analysis of the radial brightness profiles 
in terms of morphological components, reaching the following conclusions: 

• The radial light profiles show different extensions for the emission at the different wavelengths. 
The profiles from 12 to 100 jim. show a central peak, reflecting the nuclear concentration of 
warm dust due to the starburst activity. On the outer parts the profiles correspond to an 
exponential disk. At 180 jim. the slope is rather flat, and corresponds to exponential disk. 
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The decay of the light profile at 180 /xm is comparable to that seen in the optical (B band). 
An extended cold dust halo is revealed in the 180/xm map. Marginal evidence is also found 
at 100/xm. 

• We have identified three main structures: an unresolved nuclear component which corresponds 
to the central starburst, a bar, and an exponential disk. In addition we have identified a ring 
at 12 and 25/Ltm. 

• The scale length of the disk is about 3.6 kpc and is very similar in all the bands studied. In 
all filters (12, 25, 60, and 100 fim) we obtained the same asymptotic limiting values for the 
ellipticity and the position angle: e ~ 0.79 it 0.10 and PA ~ 52° it 7. These values are in 
good agreement with those derived from optical images. We conclude that the scale lengths 
of both the stars and the dust are similar. 

• The disk can be characterized by dust emission at two temperatures: a warm component 
with dust temperature of ~ 55 K and a very cold component with temperature of ~ 16 K. 
The first component results dust heated by star forming regions in the galaxy disk, while the 
very cold dust is heated by the interstellar radiation field. This work confirms the presence 
of a large amount of cold dust in the disk of NGC 253, which contains ~ 94% of the overall 
dust mass. The estimated total dust mass is 8.2 x 10^ Mq. 

• We could identify in all the light profiles a nuclear unresolved source, which was modeled 

as the PSF of the corresponding filter. The SED of this component could be modeled as 
blackbody emission at ~ 50K. The dust is heated by the intense radiation field produced in 
the nuclear starburst. 

• A far IR bar has also been detected and characterized. This bar had been previously detected 
in the near IR. The bar position angle (69.5°) and its extent are in agreement with the near-IR 
bar. The SED of the bar could be modeled with two components, a warm dust component 
at ~148 K and a cold dust component at ^23 K. 

• In the light profiles at 12 and 25 fim, we could identify a ring structure, showing at 3 kpc 
coincident with the end of the bar. The IR color profiles indicate the presence of PAHs at 
the end of the bar and coincident with the position of the ring. The detection of the ring had 
been previously claimed from both molecular gas (CO) and Ha imaging. 

• The total IR luminosity computed for NGC 253 is 4 x 10^°Lq. About half of the total 
luminosity corresponds to the central unresolved component, 31% from the disk and 17% 
comes from the bar. 

• The total star formation rate in NGC 253 is about 5.0 Mq/jx. About 70% comes from 
the nuclear starburst and the rest comes from the disk and the bar. Estimates of the star 
formation rates derived from Ha photometry are a factor ~ 10 below our estimate, which is 
consistent with previous estimates for optical extinction. 
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• The gas-to-dust mass ratio varies appreciably throughout the galaxy from a value 300 in the 
inner part to ~ 60 in the outermost part. Globally NGC 253 seems to be a dust rich galaxy, 
contrary to what was found in previous studies. 
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Fig. 1. — Contours of the ISOPHOT 180 /xm map overplotted on a B band image (Taken from the 
European Southern Observatory database). The instrumental resolution is shown at the bottom 
left corner. Contours represent luminosity in logaritmic scale and range from 3a to the maximum 
value. Contours are separated by a factor 1.8 from 1 to 230 MJy/sr. 
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Fig. 2. — Left panel: Down: Contours of the IRAS image at 12 /xm overplotted on a B band 
image (Taken from the European Southern Observatory database). Contours represent luminosity 
in logaritmic scale and range from 5a to the maximum value. Contours are separated by a factor 
1.8 from 1.3 to 340 MJy/sr. Up: IRAS image at 12//m with overlaid ellipse fits. Right panel: Up: 
IRAS image at 25 /um, with overlaid ellipse fits. Down: Contours represent luminosity in logaritmic 
scale and range from 5(7 to the maximum value. Contours are separated by a factor 1.8 from 1.5 
to 2100 MJy/sr. The instrumental resolution is shown at the bottom left corner of each panel. 
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Fig. 3. — Radial variation of the elliptical isophote parameters from the different filters versus 
radial distance from the center. The different lines representation is as follows: dotted line for 
12 |Ltm, 25 dash-dot-dashed line for 25 fim, dashed line for 60 /xm, dashed-dotted line for 100 fim, 
and continuous line for 180 /xm. Prom top to bottom we present the luminosity profiles, ellipticity 
and position angle. 
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Fig. 4. — Normalized luminosity profiles at 12 /xm (dashed-dotted line), 25 /nm (dotted line), 60 /nm 
(dashed line), 100/Lim (dash-dot-dashed line), and 180 /nm (continuous line). The luminosity profile 
in the B band is also shown (connected diamonds). The B data are taken from Pence (1980). Note 
the slow decay of the 180 /xm emission, comparable to the B emission, in contrast to the sharp 
profiles in the range 12 to 100 /xm. 
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Fig. 5. — Variation of IR colors along the major axis of elliptical isophotes. The data points 
are smoothed to the worse resolution in each color combination. Horizontal lines are the values 
expected for starbursts (SB) and discs (D) from Rowan-Robinson & Crawford (1989). 
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Fig. 6. — Surface brightness profile decomposition for all bands: 12, 25, 60, 100, and 180 fim profiles 
(from top to bottom). The unresolved nuclear component (dashed line) dominates the central 
emission except at 180 //m (see text). The bar (dash-dot-dashed line) is clearly distinguished up 
to 100 /xm. The disk (dotted line) always appears at the outermost isophotes. The ring (long 
dash-dotted line) is clearly seen only at 12 and 25 /xm. 
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Fig. 7. — Spectral energy distribution of the morphologieal components. The SED of the disk 
(dashed hne) is modeled as the emission of two blaekbodies with T ~ 55 K and T ~ 15 K. The 
SED of the nuclear component (dashed-dotted line) corresponds to a blackbody of T ~ 50 K. The 
SED of the bar (dotted line) is decomposed into two components peaking at T ~ 150 K and T 
~ 22 K. In the bottom panel we present the residuals between the sum of models and the total 
emission at each wavelength. There is a large residual at 12 ^m because there the emission is 
mainly attributed to PAH and very small grains which are in thermal equilibrium. 
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Fig. 8. — Variation along the semi-major axis of the accumulated mass for each morphological 
component. The disk is represented by a continuous line, the bar by a dashed line, and the PSF 
by a dashed-dotted line. The errors in the disk and PSF data are of the order of the size of the 
symbols used. The errors of the bar are however fairly large (^I.SIO^M©) and are not plotted to 
avoid making the graph look confusing. 
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Fig. 9. — Radial variation of gas-to-dust mass ratio. For each radius the dust mass is the sum of 
the dust masses of ah morphological components (PSF-|-disk-|-bar). 



